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Abstract 

We discuss the symmetries of quenched QCD with Wilson fermions, starting from 
its lagrangian formulation, taking into account the constraints needed for convergence 
of the ghost-quark functional integral. We construct the corresponding chiral effective 
lagrangian, including terms linear and quadratic in the lattice spacing. This allows 
us to study the phase structure of the quenched theory, and compare it to that in 
the unquenched theory. In particular we study whether there may be an Aoki phase 
(with parity and flavor spontaneously broken) or a first order transition line (with no 
symmetry breaking but meson masses proportional to the lattice spacing), which are 
the two possibilities in the unquenched theory. The presence of such phase structure, 
and the concomitant long-range correlations, has important implications for numerical 
studies using both quenched and dynamical overlap and domain-wall fermions. We 
argue that the phase structure is qualitatively the same as in the unquenched theory, 
with the choice between the two possibilities depending on the sign of a parameter in 
the low-energy effective theory. 
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1 Introduction 



Wilson fermions are one of the oldest and most important methods for formulating 
the quark sector of QCD on a lattice [1]. Not only have they (or improved versions) 
been used directly for many numerical computations, but they also play an important 
role in the more recent domain-wall [2, 3] and overlap [4] formulations of lattice QCD. 

Wilson fermions do not preserve the chiral symmetry of continuum QCD, and 
as a consequence the quark mass has to be tuned to recover chiral symmetry in the 
continuum limit [5]. Flavor and parity are exactly preserved, but, as observed by 
Aoki [6], these symmetries may be spontaneously broken by a pionic condensate for 
certain values of the bare quark mass. This leads to the existence of the so-called 
Aoki phase. In the theory with two flavors, a pionic condensate (which we can choose 
to point in the 3-direction in isospin space) breaks the SU(2) of isospin down to 
£7(1), giving rise to two Goldstone bosons at non-zero lattice spacing. Within the 
Aoki phase, the third pion has a mass proportional to a, the lattice spacing. In the 
continuum limit, the Aoki phase shrinks to zero width, the pions all become massless, 
and the pionic condensate can be rotated (by a non-singlet axial transformation) into 
the usual condensate associated with spontaneous chiral symmetry breaking. 

The existence of the Aoki phase was investigated by two of us [7] using the 
chiral effective lagrangian for pions in two-flavor unquenched lattice QCD (see also 
Ref. [8]). It was found that indeed an Aoki phase may occur near the continuum limit, 
if a certain low-energy constant in the effective lagrangian has a particular sign. For 
the other choice of sign, however, there was no spontaneous breakdown of flavor and 
parity, but rather a first order transition, in the vicinity of which the three degenerate 
pions have masses proportional to a. 

In this paper we attempt to extend the effective lagrangian investigation of Ref. [7] 
to the quenched theory. We have several motivations for doing so. First, the issue is 
theoretically interesting, and challenging, because of the presence of ghost-quarks in 
the lagrangian formulation of quenched QCD. While the peculiarities and pathologies 
of the quenched theory do show up in perturbative investigations of the effective 
theory [9, 10], they are more prominent in the non-perturbative analysis needed to 
study the phase structure. 1 

Second, the issue is of practical importance for domain-wall and overlap fermions. 
Both are built upon the (hermitian) Wilson-Dirac operator with a large negative 
quark mass. This operator is effectively quenched (even in dynamical domain-wall 
or overlap simulations) because its quark mass is not related to that of the physical 
quarks. As explained in detail in Ref. [11], the negative quark mass must be chosen so 
that one does not lie in or near any Aoki phase. The arguments of Ref. [11] also imply 
that, if there is a first order transition like that of the unquenched theory, one should 
not simulate near the phase transition line. The essential point is to avoid regions 
where the quenched pion correlators are long ranged. Thus it is clearly important to 

1 We recognize that the quenched theory is not a thermodynamic system in the usual sense, as 
the quarks do not influence the gluon fields. Nevertheless, the properties of quarks propagating on 
quenched gluon fields as a function of the bare quark mass exhibit the non-analyticities familiar 
from thermodynamics, so we think it appropriate to use the term "phase structure." 
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understand the phase structure of the quenched theory. 

Finally, numerical simulations provide some evidence for the presence of the Aoki 
phase in the quenched theory [12]. Indeed, results from quenched simulations led 
Aoki to make his original proposal of a new phase. 2 

More recent work has shown, however, that the situation in the quenched theory 
is different from, and more subtle than, that in the unquenched theory. In partic- 
ular, it was observed (see for example Ref. [14]) that there always appears to be a 
non-vanishing density of near-zero modes of the (hermitian) Wilson-Dirac operator 
in quenched QCD if the quark mass is in the super-critical region. 3 This implies a 
non- vanishing pionic condensate, and thus, through the Banks-Casher relation [15], 
would seem to lead to the conclusion that an Aoki phase fills the whole supercritical 
region, in contrast to what is expected in unquenched QCD. It was shown in Ref. [11] 
that this is not the case, however, if one defines the Aoki phase as that region of the 
phase diagram where Goldstone bosons associated with the symmetry breaking oc- 
cur. To see this, Ref. [11] considered quenched QCD with two flavors in the presence 
of a twisted quark mass, which explicitly breaks flavor and parity symmetry. It was 
argued that (in the limit of vanishing twisted quark mass) regions may exist where 
the condensate does not vanish because of the existence of a density of exponentially 
localized near- zero modes, without any of the corresponding long-range physics usu- 
ally associated with spontaneous symmetry breaking. This phenomenon is an artifact 
of the quenched approximation, and was shown to be consistent with the usual Ward- 
identity argument for the existence of Goldstone bosons. It turns out that, in the 
quenched case only, the Ward identity can be satisfied without a Goldstone pole even 
in the presence of a non-vanishing condensate, if this condensate arises because of 
localized near- zero modes [11]. Since the localization length of these near-zero modes 
is of order the lattice spacing, this phenomenon leaves no trace in the long-distance 
behavior of the quenched theory. 

The implication of the analysis of Ref. [11] is that, in the quenched theory, there 
is not a one-to-one relation between the pionic condensate in an effective field theory 
(which is sensitive only to long-distance physics) and that determined on the lattice. 
While a prediction of a non- zero pionic condensate in the effective field theory implies 
the presence of a lattice pionic condensate, the converse is not true. Thus numerical 
evidence of a condensate is not, by itself, pertinent to the question of the phase 
structure. What is pertinent, however, is the presence of a region in the phase diagram 
in which there are long distance correlation lengths, i.e. pion masses satisfying <C 
Aqcd- An effective low-energy theory can be used in any such region. A considerable 
body of numerical evidence indicates that the quenched pion mass does extrapolate to 
very small values. Based on this, we assume that there is a region where an effective 
theory can be used, and then study its properties. 

To carry out such a study we need a quenched effective lagrangian that is suited to 
non-perturbative investigations. A systematic approach to the quenched theory along 

2 In the unquenched theory, recent simulations suggest, however, that, with the Wilson gauge 
action, the scenario with a first-order phase transition applies [13]. 

3 Defined as the region where the bare quark mass satisfies — 8r < ma < 0. This is where the 
Wilson-Dirac operator can, in principle, have exact zero modes. 
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the lines of the continuum development of Ref. [16] was given in Ref. [9]. It turns 
out, however, that, while this approach is sufficient for setting up chiral perturbation 
theory (ChPT), the effective lagrangian given in [9] is not suitable for the study of 
the phase structure of the theory. 

In the approach of Ref. [9], a ghost quark is introduced for every valence quark [17]. 
The ghost quarks couple to the gluons in the same way as the valence quarks, and have 
the same mass, spin and flavor symmetries. The only difference is that ghost quarks 
have bosonic statistics. Their determinant cancels that from the valence quarks, 
thus providing a path-integral definition of quenched QCD. (For the extension to the 
case with both valence and sea quarks, see Ref. [18].) The ghost sector was dealt 
with only formally in Ref. [9], without regard to the convergence of the ghost-quark 
path integral. Since the ghosts are bosonic, this is a non-trivial issue. While the 
formal treatment of Ref. [9] is sufficient to develop quenched ChPT [19], a more 
careful treatment leads to a somewhat different symmetry structure of the quenched 
QCD lagrangian [20]. This different symmetry structure leads to a different chiral 
lagrangian, which, while equivalent to the one of Ref. [9] for ChPT, is also suitable 
for non-perturbative investigations. The construction of this lagrangian and its use 
to investigate the Aoki phase for the two-flavor theory are the central subjects of this 
paper. 

We should stress at the outset that our analysis is, in several respects, incomplete. 
At various stages we are forced to make additional assumptions not required in the 
corresponding analysis of the unquenched theory. While we think our assumptions 
are reasonable, it would clearly be preferable to avoid them. 

The plan of the paper is as follows. In Sect. 2 we define quenched QCD with iV 
flavors of Wilson fermions, paying careful attention to the convergence of the ghost- 
quark path integral. In order to do this, we formulate the theory in euclidean space, as 
is done in numerical simulations. In Sect. 3 we discuss the symmetries of this theory 
in detail, and then use these in Sect. 4 to construct the chiral effective lagrangian both 
in the continuum limit and including the leading effects of discretization. In Sect. 5 
we employ the resulting effective potential, as well as certain general properties of 
the quenched theory and the large- N c limit (N c is the number of colors), to analyze 
the two-flavor theory. We end with a summary and some concluding remarks. An 
appendix briefly reviews illustrative calculations of quenched small-volume partition 
functions. A preliminary account of this work was given in Ref. [21]. 

2 AT-flavor Quenched QCD with Wilson Fermions 

In order to define quenched QCD, we introduce, for every physical quark field q, a 
ghost-quark field q with the same quantum numbers (spin, flavor and color) as q but 
opposite statistics [17]. If there are N flavors, both q and q are A^-dimensional vectors 
in flavor space. In the continuum, the euclidean lagrangian for quenched QCD may 
then be defined as 
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C = q L Dq L + q R Dq R + q R Mq L + q L Mq R (1) 
+q* R Dq L + q[Dq R + q[Mq L + q R Mq R , 

where 

qh = Phq , qn = Pnq , 

~q~L = q p R , qr = qPi , (2) 

in accordance with standard conventions, while, in the ghost sector, 

ql = Piq, qR = Prc[, 

t = $PL, $ r = $Pr. (3) 

Our convention is P LjR — (1 ± 7s)/2. We take the mass matrices to satisfy M = M*; 
numerical simulations usually involve real diagonal mass matrices, with M = = 
M. 

In euclidean space, the Grassmann variables q R ^ and q L R are all independent, 
and the integration over them leads to the quark partition function which is just the 
determinant of the fermionic operator, 

det [D + MP L + MP R ) . (4) 

This result holds for arbitrary mass matrices M and M. For the integral over the 
bosonic ghost fields to converge, however, we must restrict M and M to be hermitian 
with all eigenvalues positive. (The euclidean Dirac operator D is anti-hermitian, and 
thus plays no role in the convergence.) Since we are taking M = M\ this means that 
M = M. The gaussian integral then leads to the ghost-quark partition function 

det - 1 (D + MP L + MP R ) = det ^ (D + M) . (5) 

This cancels the quark determinant, and we see that indeed Eq. (1) describes quenched 
QCD. We stress that q RL must be the hermitian conjugates of q Ry L for the integral 
over ghost fields to converge, while no such connection exists between q and q. These 
facts will have implications for the symmetries of quenched QCD, to be discussed in 
the next section. 

We would like to extend this continuum (and thus still formal) definition to lattice 
QCD with Wilson fermions. In the lattice action, D is replaced by the naive (nearest- 
neighbor) lattice discretization, which retains the anti-hermiticity of the continuum 
D. The mass term M is replaced by M + W, where M is a local lattice mass term 
and W the Wilson term, 

(q R Wq L )(x) = ~ ^2q R {x)r (U^{x)q L {x + fi) + U^x - fi)q L {x - fi) - 2q L (x)) . 
/* 

_ (6) 

The corresponding term with L <-> R, which we call W, contains r. In general we 

need to take r and f to be flavor matrices, so that we can treat them as spurion fields. 
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They will always be related by f = r*. In simulations, however, one usually takes 
r = f to be the identity matrix, and this is the choice we make for the remainder of 
this section. 

As is well known, in order to approach the chiral limit at non-zero lattice spacing 
one must work at negative bare quark mass. This is necessary because W, being 
positive semi-definite, makes a positive contribution to the physical quark mass. But 
once the quark mass is negative, so that M has negative eigenvalues, there are some 
configurations on which D + M + W itself can have eigenvalues with a negative or 
vanishing real part. (The case of vanishing real part corresponds to the so-called "ex- 
ceptional configurations.") It follows that we cannot simply take over the definition 
(1), with M — > M + W, to define the quenched lattice theory, since the ghost-quark 
integral would be ill-defined. The same holds true in the continuum if M has negative 
eigenvalues. 

In these cases we must proceed in a different way. We begin with the unquenched 
quark lagrangian for Wilson fermions with M = M, and we work in a basis in which 
the mass matrix M is diagonal, while adding an infinitesimal parity-odd mass term: 

N 

Here j is the flavor index, and the infinitesimal parameters 6j can have either sign 
independently for each flavor. We use primed fields in Eq. (7) because, in order to 
avoid the problem of non-positive eigenvalues of D + M + W, we perform a change 
of variables given by the following axial transformation: 



.7T 



Qj = ex P 

The quark-sector lagrangian becomes 



Qj » Qj = Qj ex P 



sg^{ej)ijl5 



(8) 



£ qU ark = ^2[qj(d + sgn{ej) t l5 {Mj + W^jqj + \ €j \ qfl 3 ] . (9) 
i=i 

The resulting lattice Dirac operator D ± i^§{M + W) is anti-hermitian, 4 and thus is 
suitable for extension to the ghost sector. Furthermore, the €j terms (which, after 
the axial transformation, look like regular mass terms) ensure convergence once we 
extend to the ghost sector. Note that for this to be true the direction of the axial 
rotation in Eq. (8) must be correlated with the sign of the infinitesimal parity-odd 
mass term in Eq. (7). 

We can now define quenched lattice QCD with Wilson fermions by adding the 
ghost sector: 

C w = Yl ( D + sgn(e J )i 75 (M j + W)) qj + ^{d + sgn(e J )i 75 (M j + wfjqj 
j'=i 

+ M^ + gk-)! • (10) 



4 This form is unfamiliar, but is, in fact, unitarily equivalent to i times the well-known hermitian 
Wilson-Dirac operator, H w = 75 [-D ± (M + W)]. 
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Because of the exceptional configurations mentioned above, quenched correlation 
functions diverge, in general, in the limit €j — >• 0, since zero modes of D + M + W are 
also, after the axial rotation, zero modes of D ± 275 (M + W) [22]. We therefore keep 
the €j non-zero, but infinitesimal. 

Several comments are in order. First, one may worry that the transformation 
of Eq. (8) could be anomalous. This is not the case. The fermionic measure on 
the lattice is rigorously invariant, and no vacuum angle 6 is generated. This is in 
accordance with the fact that in order to produce a non-zero vacuum angle from 
Wilson fermions, one has to introduce a relative phase between the Wilson mass (W) 
and the single-site mass (M) [23]. 5 

Second, the fact that the sign of €j may be chosen independently for each quark 
flavor has interesting consequences. For example, with two quenched flavors, one 
may choose e = e u = — > for the up and down flavors. This corresponds to 
(quenched) twisted-mass lattice QCD [6, 25]. It leads, after the axial transformation, 
to the lagrangian 

£ 2 w flav ° r = q(p + irsl 5 (M + wj)q + q*(D + tr 3l5 (M + W)")q + e(qq + $q) , (11) 

where r 3 is the diagonal Pauli matrix acting in flavor space. We see that, in order to 
guarantee convergence of the ghost-quark integral, a non-trivial flavor dependence is 
introduced into the quark-mass and Wilson terms. If we choose e u and of the same 
sign, however, no such flavor dependence appears. We will analyze both cases in this 
paper. 

For later use, we rewrite the lagrangian (10) in two other ways. For the sake of 
clarity we give the expressions only for the simplest choice of the €j, namely that all 
are equal and positive: €j = e > 0. It is straightforward to generalize the expressions 
to any other choice for the 6j by going back to Eq. (10). We first rewrite Cw m terms 
of the 75 projected fields of Eqs. (2,3): 

C w = q L Dq L + q R Dq R + q R Mq L + q L Wq R (12) 
+q R Dq L + q[Dq R + q[Mq L + q R Wq R , 

with M = i(M + W — ie), and e implicitly multiplied by the identity matrix in flavor 
space. Next, we group the quark and ghost-quark fields into a "super field" 

*l,r= (g'*) , *R,L = (q R , L qln) , (13) 

in terms of which 

C w = V L DV L + ^rD^r + VrMVl + VlMVr , (14) 

with 

and M = 

5 See Ref. [24] for an alternative method. 
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M =i 



M + W 




ie 





M + W -ie 



(15) 



3 Symmetries of Continuum Quenched QCD 



In this section we consider the symmetries of quenched QCD in the continuum limit. 
In that limit, the Wilson term is irrelevant once the additive renormalization of the 
quark mass has been included. Thus we can treat the quantity M + W as simply a 
quark mass matrix. Note that, because of our use of the axial transformation (8), we 
can treat the quenched theory for either sign of the quark masses. 

It is instructive to first discuss the quark and ghost sectors separately. For M — 0, 
the quark part of Cw is invariant under 

Ql,r -> V L , R q L ,R , q LtR -> q L>R V£ l R , (16) 

with Vl and Vr both elements of GL(N) (for iV flavors). Thus the chiral symmetry 
group is GL(N) L x GL(N) R . (This ignores the anomaly, to which we return below.) 
This can be maintained as a symmetry of the full quark lagrangian if M and are 
treated as independent spurion fields. Renaming the latter as M t — > M, the required 
transformations are 

M -> VrMVl 1 , M -> V L MV R L . (17) 

We take the vector subgroup to be that which remains when M and r are proportional 
to the identity matrix and all 6j are equal. This requires Vl = Vr, so that the vector 
subgroup is the diagonal subgroup GL(N). 

The symmetry group of the quark sector in the euclidean formulation is larger 
than that of the hamiltonian formulation, where the group GL(N) is reduced to the 
unitary group U(N). It is straightforward to show that the larger group (which is 
just the complexification of U(N)) does not lead to any new Ward identities (see for 
instance Ref. [19]). 

In the ghost sector, the symmetry transformations of (Jl and q R are coupled, and, 
setting M = in Eq. (12), the ghost lagrangian is invariant under [20] 

q L ^Vq L , qR^V^qR, (18) 

with V G GL(N). In other words, if c[l,r Vl,r<][l,r, the form of the ghost lagrangian 
leads to the requirement that 

V L = Vt 1 = V. (19) 
The transformations of the spurion fields M and M are 

\i ■ r 1 . .\i — r.ur . (20) 

The vector subgroup, which we define as the subgroup which leaves q L qi and 
q^R^R invariant, has VW = 1, i.e. V G U(N). If we parametrize V G GL(N) as 

V = exp(ao + a i T i) > ( 21 ) 
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with Tj the hermitian generators of SU (N) and a 0;i complex, then the vector subgroup 
has a j imaginary. "Axial" transformations are those with a and ctj real, so that 

v = vY 

Next, we discuss the symmetries of the full lagrangian, considering also transfor- 
mations of quarks into ghost quarks and vice versa. For this, it is useful to start from 
the form of the lagrangian given in Eq. (14). If \l> and \1/ were independent, then Cw 
would be invariant under 



L.R i 



(22) 



with Vl,r £ GL(N\N), as long as we treat M. and M. as independent spurion fields 
transforming as 



M -> V R MV L 



(23) 



This would imply that the (massless) quenched theory possesses the graded chiral 
symmetry group GL(N\N) L x GL(N\N) R . This does not take into account, however, 
that \& l,Ri when restricted to the ghost sector, is not independent of ^l,r, but rather, 
that 

^L,R | ghost,body = *ij,L I ghost,body • (24) 

What this equation means is the following. After a graded transformation, the ghost 
field q picks up terms containing Grassmann numbers. In general, q can be split into 
its "body" and a "soul," where the body is the usual complex field, and the soul is 
even in Grassmann fields or parameters. As discussed in more detail in Ref. [19], it 
is sufficient for the convergence of the ghost integral that Eq. (24) apply only to the 
body of the ghost fields, as already indicated in that equation. It can then be shown 
that the restriction of Eq. (19) generalizes to [19] 



Lgg I body 



Vt- 1 ! 

v Rgg I body 



(25) 



where, in N x N block form (for N flavors), 



V L qq V Lqg 



V 



Lgq ^Lgg 



(26) 



and similar for Vr. The label q refers to the quark sector, and the label g to the ghost 
sector. Note that the resulting set of transformations, 



Q' = {(Vl,Vr) e GL(N\N) L x GL(N\N) R 



V 



Lgg | body 



v Rgg I body J , 



(27) 



still forms a group. 

To complete the discussion of continuous symmetries we review the impact of 
the anomaly [9]. In the continuum limit, transformations in Q for which sdet (V L ) ^ 
sdet (Vr) are anomalous, and should be excluded. This can be accomplished by 
restricting Vl,r to lie in SL(N\N) rather than GL(N\N). There is one subtlety, 
however. The non-anomalous vector U(\) transformations do not commute with 
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general elements of SL (N\N) and so the complete symmetry group is a semi-direct 
product [26]: 

G = {(Vl,V r ) e [SL(N\N) L x SL(N\N) R ] x U(l) v | V Lg9 | body = V^| body } . (28) 

The vector subgroup (the subgroup which leaves invariant) is not affected by the 
anomaly and is given by 

H = {(V = V L = V R ) e GL(N\N) | Vj body = V^lbody} . (29) 

For further discussion of the symmetry group, see Ref. [19]. 

Finally, we consider the properties of Cw under parity. Unquenched QCD with 
Wilson fermions in the standard form is invariant under parity, but this is not obvi- 
ously true for our definition of quenched QCD. In terms of the new quark variables 
of Eq. (8), a parity transformation acts on the spinor indices as 

q(t, x) -> -x) , q(t, x) -> q(t, -0)17574 , (30) 

and the quark sector of Cw is invariant under this transformation (if the gauge field 
is transformed accordingly). The ghost sector is not, however, invariant under this 
symmetry. If we try to define parity on the ghost field as 

q(t, x) -> i7475?(^, -x) =>• <t(t, x) -> -<?(t, -x)i7 5 7 4 , (31) 

the ghost part of Cw transforms into minus itself. 

However, as long as we consider correlation functions involving only physical 
quark fields parity is not broken. Parity breaking can then only come from ghost loops, 
but these are always exactly cancelled by physical quark loops. A more formal version 
of the argument follows when we couple only the physical quarks to external sources. 
After integrating over the quark and ghost fields, the determinants cancel, and only 
the physical quark propagator couples to the sources. We conclude that parity is not 
broken in quenched QCD, with the proviso that we only consider correlation functions 
made out of physical quark fields. 

It will be useful to consider also "naive" parity: 

q(t,x) -> 7 4 <?(t, -x), q(t, x) -> q(t, -x )7 4 , (32) 
q(t, x) -> 7 4 g(t, -x) , q(t, x) ] -> q(t, -x) t 7 4 , 

Under this transformation, the kinetic terms of both the quark and ghost sectors are 
invariant, but not the mass terms, which transform into minus themselves because of 
the 75 they contain. The lagrangian (14) is, however, invariant if we transform the 
spurion fields M. and M. under naive parity as 

M^M, M^M. (33) 
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4 Quenched Chiral Effective Lagrangian 



In this section, we discuss the construction of the effective theory for the Goldstone 
particles of quenched QCD, assuming that the chiral symmetry is spontaneously bro- 
ken and that Goldstone-like excitations occur. These assumptions are based mainly 
on the results from numerical simulations. The order parameter for symmetry break- 
ing is the quark condensate, which, for positive degenerate quark masses, is aligned 
so that the vector symmetry is unbroken. 

We will need to address several issues that arise in the quenched theory that 
are not present when constructing the effective theory for unquenched QCD. The 
first is the fact that, in the quenched theory, a non- vanishing condensate does not 
necessarily lead to Goldstone-like excitations [11]. The Ward identity can be saturated 
by localized near- zero modes of the lattice Dirac operator. According to the conjecture 
of Ref. [11], such near- zero modes are present throughout the supercritical region, 
except where there are extended near-zero modes. In the latter case, the Ward identity 
is saturated by the usual Goldstone particles. The net effect of these observations 
is that the condensate in the chiral lagrangian cannot be identified with that at 
the quark level, but rather with a "long-distance" condensate defined by removing 
the contributions of short- distance near-zero modes. Since the €j — > divergences in 
quenched correlation functions are due to localized near-zero modes [11], we conjecture 
that the limit ej — > can be taken after the removal of these short-distance near- zero 
modes. 6 

The second issue concerns the presence of infra-red divergences in the chiral limit 
which make this limit singular [9, 10]. To avoid these, we assume that there is a region 
with small enough (but non-zero) quark masses for which the correlation lengths of 
Goldstone correlation functions are much larger than those of other excitations in the 
theory, so that it makes sense to use the effective lagrangian approach. 

The third issue is the need to keep the "singlet" Goldstone field, $o (defined 
precisely below), despite the fact that the corresponding symmetry is anomalous. 
This point is explained in Ref. [9, 19, 10]: $ must be kept in the effective theory 
because its correlators have long distance contributions, even though some of these 
contributions are not those of a standard single-particle pole. 

The remaining issues will be discussed as they arise in the rest of this section. In 
the first subsection we discuss the effective theory for quenched QCD in the continuum 
limit. In the second subsection we extend this to include terms of order a and a 2 , 
where a is the lattice spacing. 

4.1 Continuum Effective Lagrangian 

Here we construct the effective lagrangian for quenched QCD in the continuum limit, 
starting from the formulation of quenched QCD developed in Sect. 2. This was done 
before in Ref. [9] , but there it was naively assumed that the full chiral symmetry group 
is the graded group U(N\N) L x U(N\N) R , an assumption which we have seen above 

6 In principle, this can be done by allowing only admissible gauge fields [27]. 
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not to be entirely correct. It turns out that the effective lagrangian is equivalent to 
that of Ref. [9] if the aim is only to develop chiral perturbation theory for quenched 
QCD, but not if one wants to do non-perturbative calculations, such as needed to 
explore the Aoki phase. 

Following the standard development, we expect the Goldstone excitations to be 
described by a non-linear field E which transforms as ^l^r and spans the coset Q/H: 

E = exp($) , E -> VlEVk 1 . (34) 

Here Vl,r are elements of the symmetry group Q, Eq. (28), and H is the unbroken 
subgroup of Eq. (29). If (E) is proportional to the identity matrix, then $ is a linear 
combination of the broken generators. To construct a lagrangian invariant under Q, 
we also need E -1 — > VrE~ 1 V^ 1 , transforming as ^r^l. As already noted above, 
however, we must enlarge E to include the "super- r/" field $ = — i str log E = 
—% logsdetE (str is the supertrace, sdet the superdeterminant). This can be done 
by allowing Vl,r to be elements of the larger, anomalous, group Q', Eq. (27), while 
still only requiring the lagrangian to be invariant under Q. Since $0 is invariant under 
Q, this leads to the usual presence of an arbitrary set of functions of this singlet field 
in the effective lagrangian [16, 9]. 

With these ingredients, and treating the masses as spurions as described above, 
the 0(p 2 ) effective lagrangian becomes 

£ cff = lf 2 V l (%)sti(d^d^- 1 )-vV 2 (%)sti(M^ + ^- 1 M) (35) 

o 

+I Co $2Vb($ ) + \<x{d^ Q ) V 5 ($„) , 

with Vi($o) = 1 + 0(<&q), i = 0, 1,2,5. We have used a field redefinition to remove 
possible terms linear in $0 [16, 9]. The constant Co is of 0(A 2 ) (A = Aqcd being 
the non-perturbative scale of QCD) and, because of its relation to the rj' mass and 
topological susceptibility, is expected to be positive [9, 10]. The mass matrix M, 
Eq. (15), becomes (taking all €j positive and equal for simplicity) M. — i(m — ie) with 
m = M — m c (r) the subtracted mass matrix, 7 and Ai = Ai^ . The plus sign in the 
mass term follows from the fact that under naive parity E — > E _1 and Ai <-> Ai, cf. 
Eq. (33). Note that, while this looks like the standard mass term in ChPT, there are 
factors of ±i hidden in M. and M, related to the original singlet axial field redefinition 
of Eq. (8). 

We now consider how to parameterize the field <3>. Were it not for the restrictions 
on the ghost-ghost part of the transformations, $ would generate GL(N\N) and thus 
be a general complex graded matrix. The restriction (19) implies, however, that the 

7 m c (r) is the critical quark mass. As one approaches the continuum limit, this can be defined 
as the value for which the pions become massless. This definition does not work away from the 
continuum limit, however, since the subsequent analysis shows that the pions may have a minimal 
mass of O(a), or there may be a region of quark masses for which there are massless pions. For our 
purposes, however, it is sufficient to have a definition which determines m c to within an accuracy of 
0(aA 2 ), and one such definition is that the pion mass should be of size ah 2 . 
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ghost-ghost part of $ must have an hermitian body, and so we parameterize it as 



with lj2 hermitian N x N matrices of c-numbers, x an d X N x N matrices of 
independent Grassmann variables, and an hermitian c-number N x N matrix. Note 
that, in addition, 0i 5 2 as well as can have arbitrary souls. 

There are a number of unusual features of the parameterization (36). First and 
foremost, the presence of both 0i and 02 in Eq. (36) implies a redundancy in the 
number of fields describing the physical pions. This redundancy is not, however, 
special to the quenched theory — it occurs also in unquenched QCD. The standard 
prescription is to set 02 = 0, so that E is unitary in the quark sector, reflecting the 
fact that only the unitary subgroup are symmetries of the hamiltonian. We will follow 
this prescription here, and briefly review the arguments supporting this choice for the 
quenched theory. 8 

First, we note that the effective lagrangian depends only on the field = 0i — ?0 2 , 
and not on 0i -M02- If one is only interested in developing chiral perturbation theory 
one may just expand in (as well as 0, x an d x)- This leads to the standard 
Feynman rules in the physical meson sector. 9 However, if one wishes to consider non- 
perturbative issues such as the phase diagram of the theory, or the exact finite- volume 
partition function, the contour of integration of needs to be specified. 

We note in passing that the absence of a factor of % multiplying in (36) means 
that the propagator for ghost-ghost mesons has the sign of a physical meson propaga- 
tor (due to the minus sign from the supertrace). This differs, superficially, from the 
standard Feynman rules of quenched chiral perturbation theory in which the ghost- 
ghost mesons have wrong sign propagators [9]. There are also, however, additional 
factors of % in vertices, and it is straightforward to see that these factors can be shuf- 
fled from vertices to propagators in such a way as to reproduce the standard rules. 
This is why the use of the correct symmetry group is unnecessary when only develop- 
ing perturbation theory. Of course, in order to see the relation to standard quenched 
chiral perturbation theory one must also "undo" axial rotations of Eq. (8). 

A stronger argument for setting 2 = is that, if we choose to integrate the 
field along the contour 2 = 0, the quenched small-volume partition function is 
independent of the quark mass, as it should be [28, 29]. 10 A different contour would 
not lead to the same result. In the Appendix we give a demonstration of this in the 
one-flavor theory in the sector with zero topological charge, using a parameterization 
which makes the calculation particularly simple. 

Another feature of the parameterization (36) is the presence of the "souls" (nilpo- 
tent parts) of the fields 0i, 02 and 0. In particular, the soul of need not be hermitian, 

8 For a parallel discussion in the partially quenched case see Ref. [19]. 

9 In fact, we are assuming implicitly in this discussion that we should use the lagrangian (35) as 
is, rather than taking its real part. We are also assuming that the parameters /, v, etc. are real. 
10 Here "small volume" refers to the regime m^L « 1 « AqcdL, with L the linear size of the 




(36) 



volume. 
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although the body must be. As discussed in Ref. [19], however, the souls do not con- 
tribute to integrals, as long as the integrals are convergent. Thus, in effect, we can 
treat <f>i, 2 and as standard hermitian matrices, with the prescription discussed 
above setting 2 = 0. 

One might be concerned that if we restrict the degrees of freedom in E by setting 
2 = 0, and thus reduce the set of allowed transformations of E, we might have to 
allow additional terms to appear in the effective lagrangian. We suspect that, in 
fact, the correct approach is to use the full symmetries to determine the effective 
lagrangian, and then restrict the manifold of E. If so, this would remove the concern. 
But, in any case, the remaining transformations are sufficient to forbid additional 
terms. These transformations are unitary in the quark-quark block, hermitian in the 
ghost-ghost block, and arbitrary in the quark-ghost blocks. 

4.2 Order a Effects 

In this section, we consider how the flavor-symmetry structure of quenched QCD with 
Wilson fermions is modified at non-vanishing lattice spacing. From now on, we will 
always choose all physical and ghost quark masses equal to m. In the unquenched 
case, it was argued in Ref. [7] that the relevant continuum quark lagrangian including 
terms of order a (the lattice spacing) is given by 

£ quark , cfr = £ g iuons + q(D + i^m)q + agM 7 5 ia^F^q + 0(a 2 ) , (37) 

with b\ a coefficient which depends on the QCD coupling constant. The term linear 
in a is the Pauli term. The matrix 275 appears in both the mass and Pauli terms 
because of the field redefinition Eq. (8). The Pauli term breaks chiral symmetry in 
the same way as the mass term, and this symmetry breaking may thus be represented 
by a spurion field A which transforms just as M in Eq. (17). 

In the quenched case, analogous arguments lead to the quark effective lagrangian 

^quSS = ^gluons + *P + *75™)* + ^M 7 5 l^F^ + M + 0(a 2 ) , (38) 

where b\ can take a different value in the quenched case. We have included the 
convergence term, taking €j = e > for all flavors. Note that the effective Dirac 
operator D + 275m — bia^a^F^ is anti-hermitian, as is the corresponding operator 
in the underlying lattice theory. 

The extra 275 appears in the mass and Pauli terms because of the unusual form 
(30) that parity takes in our formulation of quenched QCD. In fact, the situation is 
somewhat subtle. Parity is broken in the ghost sector (cf. Sect. 3), and indeed, the 
ghost part of the Pauli term as well as that of the mass term in Eq. (38) breaks parity: 
like the rest of the ghost lagrangian, it changes sign under (31). However, as already 
observed at the end of Sect. 3, if we restrict ourselves to consider only correlation 
functions involving physical quark fields, parity is not broken by the underlying lattice 
theory, and correspondingly also not by >C^^ h gg. A term of the form ^fa^F^^ 
cannot appear, because it would break parity in the physical sector. Under naive 
parity (cf. Eq. (32)) the Pauli terms also switch sign, just like the mass terms. 
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It follows that, as in the unquenched theory, the Pauli term breaks chiral sym- 
metry in the same way as the mass term. We may introduce spurion fields A = ib±, 
A = —ibi = A^ transforming just like M. and M. in Eq. (23) to keep track of this 
symmetry breaking at the level of the Goldstone-meson effective lagrangian. 

With this new spurion field, we can extend the effective lagrangian of Eq. (35) 
to include 0(a) effects. The rule is simply that every appearance of M could be 
replaced by A. We are interested in vacuum structure and so will consider only the 
potential. We keep all terms of size m 2 , ma and a 2 — 0(p 4 ) in the usual chiral 
counting. Recalling that all masses are degenerate, we find 

V(S) = -icistr(S -S" 1 ) -estr(S + S- 1 ) + i c $o (39) 
+c 2 ((strE) 2 + (strE^ 1 ) 2 ) + c 3 (str E)(str 5T 1 ) + c 4 (str (E 2 ) + str (E~ 2 )) , 

with 

ci = ajmA 3 + a 2 aA 5 , (40) 
c 2 = (3 im 2 A 2 + (3 2 maA 4 + (3 3 a 2 A 6 , 
C3 = 7im 2 A 2 + 7 2 maA 4 + 73a 2 A 6 , 
c 4 = 5 x m 2 A 2 + <5 2 maA 4 + 5 3 a 2 A 6 , 

Here a iy 7, and Si are dimensionless constants of order one, and the factors of A 
are as required by dimensional analysis. The q terms in (39) can also be multiplied 
by $0 dependent potentials, but we do not make this explicit since we will argue in 
the next section that we can set $0 — in the vacuum. 

There is an apparent inconsistency in our analysis. We have kept 0(a 2 ) terms 
in the effective lagrangian, but not in the underlying quark lagrangian, £qu^keff- 
The missing terms (e.g. four-fermion operators) do not, however, break any further 
symmetries and so do not lead to additional terms in V(E). The only exception is 
a three derivative term ~ fy^^D^?, which breaks the rotation group down to its 
hypercubic subgroup, but this maps into a term with four derivatives in the effective 
chiral lagrangian, and is thus absent from V(E). 

Finally, we address the issue of the relative size of the coefficient c\ and the 
coefficients 02,3,4. Following Ref. [7], we distinguish three regimes: 

1. The quark mass is physical while a — > 0. This means m/A is a (small) constant, 
so that, for a — > 0, c 2i3i4 ~ are small compared to c\. Both lattice artifacts 
and the contributions of 0(m 2 ) terms can be ignored. In this case c\ oc m, a 
result we use several times below. In fact, the proportionality constant (ctiA 2 ) 
is very likely to be positive, 11 and we phrase subsequent discussions as though 
this is true, although our final conclusions would not change if a± were negative. 

2. The quark mass is 0(a) itself, i.e. am ~ (aA) 2 . In this case, c 2)3i4 are still small 
compared to c±, but the 0(a) term in C\ cannot be ignored. The critical value 



11 In the unquenched theory ol\ (which is related to the Gasser-Leutwyler parameter Bo) is known 
to be positive, and we do not expect quenching to change the sign of this parameter. Numerical 
evidence supports this expectation. 
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of m is shifted by an amount ~ aA 2 , and one may define m! as a subtracted 
quark mass such that the pion mass vanishes at m' — 0. The two terms in c\ 
are both of the same order. 

3. The subtracted quark mass is of order am! ~ (aA) 3 . In this case, we have that 
Ci ~ m'A 3 and 02,3,4 ~ a 2 A 6 ~ c\. In other words, for this case the higher-order 
terms in V(E) compete with the lower-order term, and we need to consider 
both to determine the vacuum structure of the theory. As already observed in 
Ref . [7] , this leads to the prediction that the width of a potential Aoki phase is 
of order a 3 (for small a). In this regime, higher-order terms beyond those shown 
in Eq. (39) are of higher order in a, and therefore need not be included in the 
analysis [7]. 

In summary, the quenched approximation has introduced several new features to 
the potential: traces have become supertraces; there are three 0(p 4 ) terms in the 
potential rather than one, due to the different group structure; there are extra terms 
involving $0; and, finally, the requirement to do an axial rotation has led to the 
leading mass-like term have the form E — E -1 rather than E + E _1 . We now turn to 
the question of how these changes influence the phase diagram. 

5 Phase Diagram 

We now address the central question of this paper, whether there can be an Aoki 
phase in the quenched theory, analogous to that predicted for the unquenched two- 
flavor theory. This might seem to be just a matter of minimizing the potential we have 
constructed, but there are a number of subtleties in the ghost sector that complicate 
the analysis. In fact, we will need certain additional properties of the quenched theory, 
which we collect before proceeding to the full analysis. 

5.1 Additional Properties of the Quenched Theory 

To study the spontaneous breaking of parity and the flavor symmetry group 7i of 
Eq. (29), we will only allow source terms which respect the graded symmetry between 
a quark and its corresponding ghost quark. The reason is, of course, that ghost quarks 
are not present in numerical simulations, but instead the quark determinant is simply 
omitted. The ghost-quark is just a field-theoretical trick to describe this procedure in 
terms of a path integral, and our choice of sources should be restricted accordingly. 
Furthermore, it is sufficient to use source terms which are diagonal in flavor. This 
is because we can choose any flavor-breaking condensate to point in the T3 direction. 
The fact that the source terms are flavor- diagonal then implies that quarks of different 
flavors are not coupled (either by the action or the sources). This in turn means that 
quenched correlation functions involving only one flavor of quark depend only on the 
mass of that quark, and not on those of other quarks, or even their presence or absence. 
It also follows that correlation functions involving ghost quarks are identical, up to a 
possible overall sign, to those for the corresponding quark, since they are composed 
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of the same propagators. These results hold for any lattice spacing and thus also in 
the continuum limit. 

A consequence of these observations is that the quenched condensate for each 
flavor is independent of the number of flavors, and that any ghost condensate is equal 
to that of the corresponding quark. This holds not only for the bare lattice conden- 
sates, but also for the physical condensates. The latter are obtained by subtracting 
divergent contributions from the bare condensate, and then multiplying by a match- 
ing factor. These steps maybe difficult to implement in practice, but what matters to 
us here is that they can be done in principle, and that the subtractions and matching 
factors do not depend on the number or properties of flavors other than that in the 
condensate itself. It follows in particular that any spontaneous symmetry breaking 
pattern does not break symmetries in Tt such that any ghost-quark condensate would 
be different from the corresponding quark condensate. 

Our next observation concerns the quenched condensate in the continuum limit. 
In this limit, one can ignore the b\ term in Eq. (38), and the quark condensate is an 
odd function of the quark mass, configuration by configuration: 

tr [1/(D + i 75 m)] = tr [75 / (D + i lb m)] = -tr [1/(D - i 75 m)] . (41) 

If the condensate has a non-vanishing non-perturbative contribution, as we are as- 
suming, then this too must be odd in m. The same holds for the corresponding 
ghost-quark condensate. Such a behavior is indeed seen in numerical results for 
the quenched condensate with staggered and overlap fermions, which is odd under 
m — > —m. At the level of the effective chiral lagrangian, working in the continuum 
limit means that c\ oc m and 02,3,4 oc m 2 in Eq. (39), and the result (41) implies 
that Egg (the component of the condensate for quark q), and the corresponding ghost 
condensate should, in the continuum limit, both change sign when c\ does. This 
result also follows from symmetries in the following somewhat indirect way. First, for 
the condensate in the physical sector, Egg, this can be derived using the appropriate 
axial symmetry in order to flip the sign of m. Then, because of the arguments given 
above, the corresponding ghost condensate, has to follow suit. We note that 
there is no "direct" way of changing the sign of m in the ghost sector, as follows 
from Eq. (20). We note also that this result does not hold in the unquenched theory 
because the measure depends on the quark mass. 

In light of the previous discussion, one might wonder why we need to consider 
more than one flavor in the quenched theory. The reason we need two flavors is to 
allow the calculation of flavor non-singlet pion propagators. These have only quark- 
connected contributions and are the quantities usually calculated in simulations. It 
is the masses of these pions which are observed to extrapolate to zero at non-zero 
lattice spacing — the phenomenon which the Aoki phase was introduced to explain. 
With one flavor alone one can only consider flavor singlet pions, which have quark- 
disconnected contributions as well. 

We also need to recall some properties of the quenched theory in the limit that 
N c , the number of colors, becomes large. 12 In the continuum effective lagrangian, 

12 One subtlety of this limit is that we must take N c — > 00 before m — > 0, so that the quenched 
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Eq. (35), standard arguments show that the constants f 2 and v are proportional to 
N c , while c and a are 0(1). The latter two are suppressed because they multiply 
terms corresponding to disconnected diagrams at the quark level. Similarly, in the 
potential for non-zero lattice spacing, Eq. (39), the coefficients of the single-supertrace 
terms, c\ (which is proportional to vm) and C4, are 0(N C ), while the coefficients of 
terms with two supertraces, c , c 2 and c 3 , are 0(1). Thus in the large-iV c limit we 
can set c , c 2 and c 3 to zero. 

The advantage of this limit is, of course, that the quark sector becomes physical, 
since the quark determinant is suppressed by 1/N C . In other words, there is no 
distinction between quenched and unquenched theories in this limit. Thus the results 
of our analysis of the quenched theory must, in the quark sector, coincide with those 
of the unquenched two-flavor theory in this limit. The analysis of the Aoki phase 
in the unquenched theory in the large-iV c limit is a simple generalization of that at 
finite N c [7] . In particular, the competition between two contributions which leads to 
the non-trivial phase structure remains in the large- N c limit (since these are the c\ 
and C4 terms). The only significant change is that the flavor group SU(2) is enlarged 
to U(2). This larger symmetry group may then be broken down to U(l) x U(l) (as 
opposed to a single U(l)) in the Aoki phase, so there are still two exact Goldstone 
bosons. Further details will be discussed below. 

The particular utility of the large- N c limit is that it allows us to check our method 
of analysis of the ghost sector. We have argued above that, with sources appropriately 
chosen, the quark and ghost condensates must be the same. Since the analysis in the 
quark sector becomes unquenched as iV c ^ 00, this means that we know the result 
which we must obtain in the ghost sector in this limit. 



5.2 Phase Diagram as N c — >■ 00 

We start from the effective potential V(£) of Eq. (39) for two flavors, u and d, and 
substitute the form for $, Eq. (36). In fact, we can simplify this form as follows: 
first, by setting x — X — 0, since these are the solutions of the classical equations of 
motion for these fields; second, by setting <fi 2 = 0, as discussed in the previous section; 
and, third, by keeping only flavor-diagonal components of the fields since we force 
any flavor-breaking to lie in the T3 direction. Thus we can investigate the vacuum 
structure using an expectation value 

(£) = £ = diag(e^", e i(t>d , e^,e* d ) , (42) 

where —71 < 4> u ^ < n are real phases, while 4> u ^ are real variables to be integrated 
along the entire real axis. 

In this subsection we consider the large- N c limit, and thus set Co = c 2 = c 3 = 0. 
The potential is then 



V = 2£ 

j=u,d 



sgn(ej)ci sin(0j) — \ej\ cos(0j) + C4cos(20j) (43) 



artifacts proportional to ln(m) /N c vanish. Since we are implicitly working at small but non- vanishing 
quark masses, so as to avoid these artifacts, we are able to take the large- N c limit without problems. 
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+i sgn(ej)ci sinh(0j) + \ej\ cosh(0j) — C4 cosh(20j) . 

Here we have reintroduced the dependence on the e^, which follows from the fact that 
the lattice mass and Wilson terms in Eq. (10) are proportional to sgn(ej). The c 4 
term does not depend on the sign of ej since it is quadratic in symmetry breaking. 
The factorization into four parts, dependent respectively on <p u , (pa, 4> u and <pd, is the 
result of taking N c — > 00. This is how, in this context, the quark and ghost sectors 
decouple. 

In the quark sectors the analysis is essentially a recapitulation of that of Ref. [7] 
for the unquenched two-flavor theory, except that here it is done for each flavor 
separately. The notation is also different: our C4 was called — c 2 /4 in Ref. [7], and 
here we have done an axial transformation which complicates the interpretation of 
the condensate. What remains unchanged is the region of parameters of interest: c\ 
is proportional to the 0(a) shifted quark mass m! while C4 ~ a 2 is a constant. A 
useful ratio is rj = ci/Ac^. We want to determine what happens to the condensate as 
i] ranges from values smaller than —1 to larger than 1, in the limit e — > 0. We note 
that the term proportional to \ej\ can be dropped for most of the analysis, since we 
take €j — > at the end, except where it is needed to distinguish between otherwise 
degenerate minima. 

It is useful, first, to consider the continuum theory in which C4 = 0, so that 
the potential in each quark sector is V q = 2 sgn(e ? )ci sin(0 g ), where we now use an 
index q instead of j to indicate that we are considering the quark sector of a given 
flavor q. For positive c\ (corresponding to positive quark mass in the original basis) 
the minimum is at 4> q = — sgn(e g )7r/2, so that T, qq = e % ^ q = — sgn(e g )i. For negative 
Ci, the minimum switches to (p q = sgn(e g )7r/2, so that Tl qq = sgn(e q )i. In short, the 
minimum occurs at Y> qq = —sgn(e q Ci)i. This reproduces the standard discontinuous 
dependence of the direction of the condensate as a function of the direction of the 
mass and source terms. 

The unusual factor of — i in the condensate is a reflection of the axial transforma- 
tion Eq. (8) is needed to define the quenched theory. If we rotate back to the usual 
continuum basis, then the condensate becomes E OT = sgn(ci). In other words the 
non-standard values of <p q "undo" the axial transformation. To see this explicitly, we 
can write the potential in terms of shifted fields (p q = — sgn(e 9 Ci) n/2 + <p': 

V q = -2| Cl |cos0'-2c 4 cos20' = -2(| Cl |+c 4 ) + (|ci|+4c 4 )0' 2 + O(0' 4 ) , (44) 

where we have restored the C4 term in anticipation of the discussion below. This is 
the potential we would have obtained directly if we had not had to use the trick with 
the axial transformation of Eq. (8), and instead expanded about S gg = ±1. 

Now we consider non-vanishing c 4 . If c 4 > 0, the analysis with c 4 = goes 
through unchanged: the c 4 cos(20 g ) term has equal minima at <p q = =F7r/2, so it is the 
sign of t q c\ which determines the direction of the condensate. There is thus a first 
order transition when c\ changes sign. 

If C4 < 0, however, it contributes negative curvature to V q , which can desta- 
bilize the minima. Looking at the quadratic terms in Eq. (44), we see that this 
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happens when 4c 4 < — |ci|, i.e. \r}\ < 1. The minimum of the potential shifts to 
4> q = sgn(e g ) arcsin(^), signaling that we are in an Aoki phase. There are two choices 
of branch for the inverse sine, one interpolating between — tt/2 and +7r/2, the other 
between — n/2 and — 3tt/2. In both cases the condensate "swings" from —% to +i 
around the unit circle in the complex plane, when e q c\ changes from positive to neg- 
ative values. The choice of branch decides which way the condensate swings, i.e. 
whether it passes through +1 or —1 when r\ = C\ = 0. Which choice is appropriate 
is determined by the \ej\ term in Eq. (43). This term is negative if — 7r/2 < <p q < n/2 
and positive if — 3n/2 < (fi q < —n/2, and thus selects the branch in which T, qq swings 
through +1. This is true irrespective of the sign of €j. 

The symmetry that is broken by the condensate depends on the details of the 
theory. If we were considering one flavor alone, then, undoing the axial transformation 
(8), we would find that (qi^q) oc sgn(e q ) a/1 — rf ^ for \r]\ < 1. Thus parity is 
broken in the Aoki phase. The sign of the condensate is determined by the sign of 
e q = €j in Eq. (7). There are no exact Goldstone bosons in this phase as no continuous 
lattice symmetry is broken. 

Next, we consider the two-flavor theory, with degenerate masses. The pattern of 
symmetry breaking depends on the signs we choose for e u and in Eq. (7). If we 
choose e u and both of the same sign, so that Eq. (7) represents a flavor singlet source 
proportional to ^2 q qij5q, then the condensates will satisfy (ui^u) = (dij^d) ^ 0, 
and parity, but not flavor, will be broken (so again there are no Goldstone bosons). 
We checked explicitly, from Eq. (39) with c = c 2 = c 3 = 0, that in this case indeed 
there are no Goldstone bosons (all pions have a mass of order a) . If, on the other hand, 
we use a flavor-breaking source term proportional to ui^^u — di^^d {i.e. €d = — e u ), 
the mass term in the lattice quark lagrangian will be proportional to r 3 (c/. Eq. (11)), 
which translates into the C\ terms in Eq. (43) for the up and down sectors having 
opposite signs. The condensates will then satisfy (uij^u) = —(di^d) ^ 0, and flavor 
and parity will be broken. This is the most familiar example of the Aoki phase, which 
has two Goldstone bosons 71^, and which carries over to finite N c in the unquenched 
theory (see Sect. 5.3 below). In this case, the calculation of the meson spectrum 
from the chiral lagrangian is virtually identical to that in Ref. [7], and confirms the 
existence of two massless Goldstone bosons. We stress, however, that, as far as the 
condensate is concerned, these theories at infinite N c differ only kinematically — the 
underlying phenomenon is identical in all cases. 

We now turn to the ghost sector, in which the potential is 



i c[ sinh(0) + |e| cosh(0) — c 4 cosh(20) . (45) 



We have dropped the index j because, in the remainder of this section, we are only 
interested in demonstrating that, for each flavor, the ghost sector follows the quark 
sector. We have introduced the useful variable 13 c[ = sgn(ej)ci. As in the quark 
sector, we can ignore the |e| term except when it breaks degeneracies. 

13 We could have performed the analysis in the quark sector in terms of c[ as well; however, there 
we wanted to emphasize the explicit flavor dependence of the condensate in the Aoki phase through 
the signs of the e q . 
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The potential (45) is complex, and the usual approach of minimizing it to find the 
vacuum structure breaks down. However, what one is really doing if one minimizes 
the potential in the context of euclidean field theory is to calculate the leading term 
in a saddle-point expansion of the path integral. Stating the problem this way, it 
is clear that this can also be done when the integrand is complex instead of real. 
The prescription is to treat as complex, and deform the contour of integration so 
as to pass through a saddle-point. The vacuum expectation value of the field is, at 
leading order, the value at the position of the saddle. If there are multiple saddles, 
the appropriate one to use is determined by the boundary conditions on the contours, 
and other considerations, as will be discussed below. 

Another issue we must face is the convergence of integrals as — > ±oo. For 
instance, if C4 > 0, the last term in V g is unbounded from below in these limits. In 
fact, we only know the form of the potential for small fields, |0| < 1. The ratio of 
successive terms in the chiral expansion grows exponentially when |0| exceeds unity 
(as can be seen from the relative magnitude of the C4 and c\ terms) rapidly overcoming 
any suppression in their coefficients. Thus, once |0| ~ 1 the chiral expansion breaks 
down, and we do not know its behavior at large fields. In light of this, we simply 
assume that we can treat the integral over as convergent at infinity along the real 
axis, and search for saddles in the region |0| <1 through which the contour can be 
deformed in such a way that it can be evaluated using steepest descent. Note that the 
same issues do not arise in the (f> q integrals, because there is no relative exponential 
growth in the terms, and because the integral runs over only a short segment (of 
length 2tt) of the real axis. 

As for the quark sector, we warm up by considering the phase structure for C4 = 0, 
so that V g = 2ic[ sinh(0). The saddle point equation is 

cosh(0) = cosh(a;) cos(y) + % sinh(:r) sin(y) = , (46) 

where we have written = x+iy, with x and y real. The solutions to this equation are 
a periodically repeating sequence along the imaginary axis, with the two nearest to 
the origin being at 4>a — —iir/2 and 0# = +in/2. These correspond to the condensate 
taking the values S 99 = exp(0A,s) = ^fi, i-e. the same values as arose in the quark 
sector. The other solutions lead to the same two values of the condensate and are 
thus not physically distinct. 

Expanding about the saddles, we find 

V g (0 = ^i 1 - + 0') = ±2c[ cosh(0') = ±2q 

The criteria we use to choose the saddle are based on the direction in which $tV g rises 
most steeply (which is the direction the contour should follow to have the steepest 
descent of | exp(— V g )\ away from the saddle), as well as the value of the potential at 
the saddle. The contour should pass through the saddle and then be able to join the 
real axis for large fields, without passing near other saddles. If there is more than 
one saddle point with an appropriate contour, we want to maximize the value of $lV g , 
so that the contribution to the saddle-point integral is minimized, and thus a better 



1/2 



1 + 



(47) 
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approximation to the (absolute) value of the full integral is obtained. 14 All criteria 
agree in this case: if c[ > then V g is larger at (f>A, and has steepest descent in the 
direction of <p' real, so that the contour can be easily deformed so as to end up on 
the real axis for large . The other saddle has its direction of steepest descent along 
the imaginary axis, with the contour heading directly for saddle <pA and its periodic 
reflection. If c[ < 0, the role of the two saddles is reversed, and we must choose <f) B . 

In this way, the ghost condensate is predicted to flip from —i to +i as c[ passes 
from positive to negative values. This is exactly the same dependence as for the quark 
condensate, H qq , and is in agreement with the general arguments given in Sect. 5.1. 
That this check works out gives us confidence in our method. We note for future 
reference that the total potential, V q + V g , vanishes at the saddles that have been 
chosen. 

Now we add back in a non-zero C4, in which case the saddle point equation 
becomes 

ic[ cosh(0) - 2c 4 sinh(20) = . (48) 
The solutions are 4>a,b, defined above, and an additional saddle or saddles satisfying 

0c : sinh(0 c ) = </(4c 4 ) = irf . (49) 

The solutions of this equation (apart from periodic repetitions in the imaginary di- 
rection) are 

7/ < — 1 : <pc± = —in/2 ± arcosh(— rf) = <\>a ± arcosh(— rf) (50) 
— 1 < 7/ < +1 : 4>c = iarcsinr/ 

7/ > 1 : 4>c± = 27r/2 ± arcoshr/ = 4>b ± arcoshr/ . 

Expanding V g around these saddles we find 

V g {h + &) = 2(c , 1 + C4)+4c 4 (V + l)0 /2 + O(0 /3 ), (51) 
V g (4> B + ft) = 2(-c / 1 + C4)+4c4(-V + l)0 ,2 + O(0 ,3 ), 
V g {4>c + ft) = -2c4(l + 27/ 2 )+4c 4 (V 2 -l)0 /2 + O(0' 3 ). 

As above, our approach will be to find the saddle for which V g increases for real 
displacements (which we will refer to as "positive curvature"), since the original 
integral over (f> is along the real axis, and which has the largest value of the potential 
at the saddle. We run through the choices of parameters in turn. 

We begin by considering C4 > 0. 

> 4c 4 > 0, or 7]' > 1: Saddles A and C± (arrayed either side of B) have positive 
curvature, while B has negative curvature. There are two choices of contour: that 
passing through A alone, and that running through (fic-, <t>B and then 4>c+- We 
choose the former as the potential is higher at A than at C±. Furthermore, saddles 
C± are outside the range of applicability of our truncated effective potential except 
for T]' ~ 1. Thus we find the same saddle as for C4 = 0. 

14 For a lucid description of the saddle-point method, see e.g. Ref. [30]. 
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4c 4 > c' x > 0, or < rj' < 1: Both saddles A and B have positive curvature, while 
C lies on the imaginary axis between A and B and has negative curvature. Thus a 
contour passing through either A or B is possible. We choose that through A since 
it has the higher potential. 

4c 4 > > c' x > — |4c4 1 , or —1 < rf < 0: The saddles are similar to the previous 
case, except that saddle B now has the highest potential, and so we pick the contour 
passing through B rather than A. In other words, when c[ passes through zero the 
saddle switches from A to B, just as in the case c 4 = 0. 

4c 4 > > — 4c 4 > c'j, or rf' < — 1: The situation is the reflection in the real axis of 
that for c[ > 4c 4 > 0. Saddles B and C± have positive curvature, and we choose the 
former since it has the higher potential. 

Summary for c 4 > 0: The ghost condensate is unchanged from the analysis when 
c 4 = 0, and satisfies <p = i<p for all c[. Thus H qq = as expected. In both quark 
and ghost sectors there is a first-order phase transition as c[ passes through zero, just 
as in the continuum limit. There is no Aoki phase in either sector. The spectrum of 
fluctuations about both quark and ghost vacua is the same for all c[. 

We now turn to the case c 4 < 0. 

c'j > 4|c 4 | > > 4c 4 , or rj' < —1: Changing the sign of c 4 changes the sign of the 
curvature about all saddles. Thus, in this parameter range, only saddle A has positive 
curvature. The contour must pass through A, and then move back to the real axis. 
It can do so by passing through C± and then moving off in an imaginary direction, 
although this is not necessary. In all cases the integral is dominated by the value at 
saddle A. 

4|c 4 | > |cj| ; > c 4 , or \r]'\ < 1: The contour must pass through saddle C, since 
this is the only saddle with positive curvature. This gives an expectation value = 
iarcsin?/ = i sgn(e.,) arcsin??, which is consistent with the equality = i(f>, since we 
found above that (fi q = sgn(ej) arcsinr/. As for the quark sector, however, there is a 
choice of branch of the inverse sine function. The appropriate choice is determined 
by the 2|e,'|cosh0| term in V g : one should pick the branch which maximizes this 
term. It is easy to see that the resulting branch is that which satisfies the equality 
<f> — i<f> for either sign of €j, i.e. that the ghost condensate is always equal to the quark 
condensate. 

> 4c 4 > c'j, or rf > 1: The contour passes through B, which is the only saddle 
with positive curvature. 

Summary for c 4 < 0: As for c 4 > 0, these results are consistent with those from 
the quark sector and give S 99 = In particular, we find an Aoki phase in both 
quark and ghost sectors. 

In summary, we see how the saddle-point analysis effectively restores the quark- 
ghost symmetry which had been broken by the need for convergence of the ghost 
integral. 
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5.3 Phase Diagram for Finite N c 



We now extend our analysis of the quenched phase diagram to finite N c . Finite- N c 
corrections will change the values of C\ and C4, but, more importantly, turn on cq, c 2 
and c 3 . The latter terms couple the quark and ghost sectors, as well as the different 
flavors. To make progress here we will need to use the assumption, discussed in 
Sect. 5.1, that the condensates in the quark and ghost sectors have to be equal, as 
well as other general considerations. 

We note first that it is straightforward to show that the solutions we found in 
Sect. 5.2 are still solutions of the full saddle-point equations which follow from Eqs. 
(39) and (42). The full saddle-point equations in the finite- iV c case are transcendental 
equations, due to the Co term in the effective potential V(E), and we do not know 
whether any solutions to these equations other than those discussed in the previous 
section exist. However, in the previous section we found that the ghost and quark 
condensates were always equal in the N c — > 00 limit, and, in accordance with the 
argument at the quark level in Sect. 5.1, we will insist that this continues to hold 
for finite N c . This means that we require <p q = i<p q for each flavor q. Remarkably, 
upon imposing this constraint we find that the terms proportional to c , c 2 and c 3 
drop out of the saddle-point equations, and the solutions are those found in Sect. 5.2. 
Note that a similar argument applies to the terms of order $q and beyond 15 in the 
potentials Vo^^o) in Eq. (35). 

We conclude that reducing N c from infinity does not change the saddles we need 
to consider. This is the good news. The bad news is that, since we are now considering 
the coupled (f> q , (f> q system, we must use the full potential when distinguishing between 
the saddles, and, as noted above, this potential vanishes for all of them. 

To make progress, we consider in general the role of the co$q term in the quenched 
chiral potential (39). A signature result of quenched chiral perturbation theory is that, 
at tree level, this term does not shift the position of the pole in the flavor-singlet 
propagator but rather gives rise to a double-pole contribution at the same position. 
This double-pole term represents the presence of quark-disconnected contributions 
(the so-called "double hairpins") in this channel [9, 10]. 16 Furthermore, the cq term 
does not affect the value of the condensate at tree level, but rather only at one-loop 
through long-distance effects. This is very different from the effect of such a term in 
the unquenched theory. There it shifts the position of the rj' pole, and changes the 
value of the condensate at leading order (since 7/ loop effects are short distance they 
lead to an 0(1) correction to the condensate proportional to v x M 2 , / (4tt f w ) 2 ) . 

We conjecture that these perturbative results hold also non-perturbatively, and 
that the c term does not influence the value of the condensate. A similar discussion 
applies to c 2 and c 3 terms in the potential, which both have two supertraces, and 
at leading order in a field expansion are proportional to $q. Given this conjecture, 
the analysis for finite N c collapses to that for infinite N c , and we conclude that the 
possible phase structures in the quenched theory at finite N c are the same as those 
in the unquenched theory at infinite N c . 

15 These terms all vanish in the large- N c limit. 

16 We expect this to be true inside the Aoki phase as well. 
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We illustrate the argument for the irrelevance of c , c 2 and c 3 by studying the 
small fluctuations around one of the saddles. For simplicity assume that c[ > and 
that c[ + 4c4 > as well, so that we are in the phase with i<p q = <p q = — in/ 2 (this 
is saddle point A in the ghost sector). Expanding the effective potential around this 
saddle to quadratic order gives 

V(E) = ±(<j> q , 4> q ) [2(ci + 4c4)I + (c + 4c 2 + 2c 3 )x] (^j , (52) 

in which I is the unit matrix and X is a matrix filled with l's in the quark sector (upper 
left-hand block), i's in the mixed quark-ghost sector, and —l's in the ghost sector 
(lower right-hand block). Note that the coefficient of I reproduces the coefficients of 
the quadratic terms in the first line of Eq. (51), as well as Eq. (44) for the case at 
hand. 

First, ignoring X, we see that the potential is stable around our solution if 
c[ + 4c4 > 0, with an instability developing when c[ + 4c4 < 0, signaling sponta- 
neous symmetry breaking. This confirms what we found in Sect. 5.2. As long as 
c[ + 4c 4 > 0, <p q and <p q mesons propagate with masses proportional to c[ + 4c 4 . 

Now if we include X, the (zero-momentum) propagator can be determined from 
the inverse of the matrix in Eq. (52). Because X 2 = 0, this inverse is equal to 

1 I- C ° + 4C2 + i C3 X. (53) 



2(ci + 4c 4 ) 4(ci + 4c 4 )< 

This shows that the parameters c 0i 2,3 do not affect the meson masses, but instead 
determine the residue of the quenched double pole, which originates in the so-called 
"double-hairpin" diagram. While this double pole exhibits a sickness of the quenched 
theory, it does not affect the meson masses. Spontaneous symmetry breaking only 
occurs when the squares of these masses turn negative. 

This brings us to our final point of this section, which is the effect of the flavor- 
singlet pseudoscalar in the unquenched theory. For simplicity we discuss only the two- 
flavor theory. In the unquenched theory, the flavor- singlet mass term is (1/2)co$q = 
(l/2)co(0 u + 4>d) 2 - In the limit of vanishing e U) d this term raises the energy of 
the solution for which <p u = <pd- Thus, in the unquenched theory, one finds that 
(ui^^u) = —(di^d) 7^ inside the Aoki phase, and flavor is always broken along 
with parity. As we have seen, this mechanism for picking the vacuum is not present 
in the quenched theory. In the quenched theory, for both infinite and finite N c , the 
pattern of symmetry breaking depends on the relative sign of e u and e^. 



6 Conclusion 

The aim of this paper is to see whether the method of Ref. [7], based on the effective 
theory describing the Goldstone-boson physics of full QCD with two flavors of Wilson 
fermions, can be extended to the quenched theory. Our conclusion is that it can, 
although to push through the analysis we have had to make a number of assumptions 
not needed in the unquenched theory. We find the same two possibilities for the phase 
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structure as in the unquenched case: depending on the sign of a parameter of the 
effective theory (c 4 ), there is either an Aoki phase or a first-order phase transition. 

The nature of the quenched Aoki phase, if there is one, can however differ from 
that in the unquenched theory. In particular, the form of the condensate in the 
two-flavor quenched theory depends on the source employed to probe spontaneous 
symmetry breaking. One can have a phase in which both parity and flavor symmetry 
are broken, as in Aoki's original scenario, but one can also have only breaking of parity, 
with no breaking of flavor, and thus no Goldstone bosons. Since in the quenched 
theory quark correlators are only probing the theory, while not being part of the 
dynamics (there are no sea quarks), this is a kinematical effect. For more detail, 
we refer to Sect. 5.2. Unquenched QCD at N c = oo shares this property, but there 
the degeneracy between the two possibilities is lifted at finite N c , in accordance with 
Ref. [7]. 

The extension of the method of Ref. [7] to the quenched theory turned out to 
be non-trivial. The formal definition of quenched QCD of Ref. [9] does not lead 
to a convergent path integral in the ghost sector. While expanding the lagrangian 
of Ref. [9] around S = 1 leads to the correct version of quenched ChPT, it breaks 
down non-perturbatively. We gave a non-perturbatively valid path-integral defini- 
tion of quenched QCD, and analyzed its symmetries. This has been discussed before 
in the continuum in Ref. [20], but here we extended this to lattice QCD with Wil- 
son fermions, thus providing a fully-regulated, convergent path-integral definition of 
quenched QCD. We then used it to construct an effective potential to next-to-leading 
order, including terms up to order a 2 . Employing this effective action, we argued that 
the analysis of the phase structure of unquenched QCD with two Wilson fermions of 
Ref. [7] carries over to the quenched case. As a by-product we saw how standard 
quenched chiral perturbation theory is recovered. 

As discussed in some more detail in Sect. 5.2, for certain values of the coefficients, 
the effective potential can be unbounded from below for large values of the fields. 
While the vacuum structure we found satisfies local stability, this would seem to cast 
a doubt on our conclusions. We believe that this problem is spurious. As we argued, 
the effective theory is only valid below the typical hadronic scale, and should thus 
not be trusted for large values of the fields. In fact, we expect that it is possible to 
construct a different effective potential with the same vacuum structure and the same 
perturbative expansion, but which is bounded for large values of the fields [31]. 

Finally, we remark that the entire analysis was carried out in euclidean space. 
This is the relevant setting: quenching is only employed in the euclidean version of 
lattice QCD. It is interesting to note that, while the euclidean ghost action is invariant 
under SO (4) rotations, the Minkowski version would not be invariant under Lorentz 
transformations, because in the ghost sector the option of identifying q with g^7° 
does not exist. Like with parity, this would not affect quenched correlation functions 
with only physical (and no ghost) quark fields on the external legs. However, we 
doubt that it will be possible even in principle to continue the quenched theory to 
Minkowski space. 
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Appendix 

First, we calculate the small- volume partition function for the one-flavor quenched 
theory in the sector with topological charge zero. The latter condition implies that 
the constant Co = 0. We include this appendix only for pedagogical reasons; for other 
cases as well as references, see for example Refs. [28, 29]. In a small volume (often 
referred to as the "e-regime"), the dominant contribution comes from the constant 
modes. The partition function we wish to calculate is 



Z(m,J) = J to exp | — v V str 



m+ J 
m 



(£ + £■ 



(54) 



\ AB B AB 2 

where V is the four-dimensional volume, and the fields are constant. A and B rep- 
resent c-number degrees of freedom, while C and D represent Grassmann degrees of 
freedom. The integration measure 2mu = —dA dB dC dD is obtained starting from 
the Haar measure on GL(1|1), and thus, following our arguments in Sect. (4.1), we 
restrict the bosonic part of the integration region to R — S± x R + [28]. In other 
words, A lies on the unit circle S 1 and B lies on the positive real axis R + . We took 
the masses in the physical and ghost sectors m + J and m different; the quenched 
theory is obtained by setting J = 0. 

A convenient parameterization for the non-linear field S equivalent to that used 
in the text is given by choosing 

A = , B = e*' . (55) 

The relation between this parameterization and that of Eq. (36) can be worked out 
order by order by an expansion in the fields. (This expansion is finite, because of 
the Grassmann nature of the fields Xi X, C and D.) The advantage of our new 
parameterization is that the jacobian which appears when we actually calculate the 
integral is simpler. 

The jacobian for the transition from A, B to 0', 0' is equal to iAB = ie 1 ^ . 
Expanding in C and D and defining r = 2vVm, s = 2vVJ, we find that 

Z(m,J) = -— d(f)' dip 1 I dCdD(iAB) (56) 



2tu 
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exp + S )(A + 1A4) - i r(B + 1/fl)} {l + ^ (l + f + l) } 

— J d(f)' J d(j)' |(r + s) cos0' + r cosh0'| 

x exp | (r + s) cos <p' — r cosh </>' j 
= (r + s)/i(r + s)Xo(r) + r/ (r + s)Xi(r) . 

To obtain the quenched result we set s = and find that Z(m, 0) = 1, irrespective 
of the value of m. 17 This would not have been true if we had chosen A to lie on R + , 
which corresponds to choosing 0i = instead of 02 = in Eq. (36). 

In the case with more than one flavor, there is more than one way to define the 
vector- like subgroup Ti of Eq. (29). For instance, for the two-flavor case, if we have a 
mass matrix proportional to r 3 , it is natural to define the vector-like subgroup in the 
ghost sector by requiring that V is r 3 -hermitian: 

V*t 3 V = r 3 , (57) 

and likewise, that Vl = t 3 Vrt 3 in the quark sector. The relevant group integral for the 
standard quenched two-flavor case, analogous to the one-flavor case reviewed above, 
was done in Ref. [32], and here we discuss only how things change if one chooses 
vector-like transformations in the ghost sector to be T3-hermitian. 

First, let us consider the parameterization of the bosonic ghost block of S, i.e. B 
of Eq. (55) above for the two-flavor case. In Ref. [32], the following parametrization 
was chosen (showing only bosonic parameters, i.e. setting Grassmann parameters 
equal to zero): 18 

Sl+S2 /cosh$ e Sl ~ S2 isinh* e iu \ , ■ 

\ i sinh $ e- ia cosh $ e -^- s ^ J ' [ } 

This can be written as 

_g _ g(0O+T"303)/2 e T 1 <t> 1 +T 2 4>2 g(0O+T"303)/2 {59) 

with 

(60) 



Sl + s 2 = 


00 , 


s l — $2 — 


03 , 


$ = 


V0? + 0i , 




01 


sin a = 


, COSCT 







02 



17 Note that the point r = s = is singular. 

18 We use denote the parameter <f> of Ref. [32] by $ here in order to avoid confusion. 
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This parametrizes the coset GL(2)/U{2) for the standard choice of H. For our new 
definition of H, the relevant coset is instead parametrized by 



£> _ e (^0+r 3 ^ 3 )/2 e iri0i+jr 2 2 e (<fo+T3<fe)/2 (g;Q 

where now the matrix in the middle is unitary instead of hermitian. Correspondingly 
we change the parametrization of the bosonic valence block, A of Eq. (55), from 

A = e Wo+ iT 3^3)/2 e iT±4>l+iT 2 4>2 e («/>O+*T-303)/2 ^g2) 

j^+i^ (cosO e^ 1 -^ 2 sinfle^ 

-sin^e-v costfe^ 1 "^ 



in which in terms of the variables of Ref. [32] 



^1+^2 = 0o , (63) 

e 



>3 



01 02 

smp = -j , cosp = — , 

to 

A T . = e (^o+* T 3^3)/2 e T 1 4> 1 +T 2 (j)2 e (i0o+*T303)/2 j-g^-j 

In other words, for those generators in which we change the fields from non-compact 
to compact in the ghost sector (0i,2), we change the fields from compact to non- 
compact in the quark sector (to be compared with the standard choice 02 = in 
Eq. (36)). 

In this parametrization the measures for the matrices W\p (parametrized by 0o ; 3 
and 0o,3) and the matrices w and w (parametrized by 0i >2 and 0i )2 ) of Ref. [32] factor- 
ize. Following section (3.1.2) of Ref . [32], our change in parametrization corresponds 
to an interchange of the variables $ and 9, which leaves the measure (cf. Eq. (3.51) of 
Ref. [32]) invariant. We thus conclude that if the partition function for the two-flavor 
case is independent of the quark mass (as we showed it to be for the one-flavor case 
above), it is also independent of the quark mass for our new parametrization, if we 
follow the prescription given above. 
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